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ABSTRACT 



We present observations of the main-belt comet 259P/Garradd, previously 
known as P/2008 Rl (Garradd), obtained in 2011 and 2012 using the Gemini 
North Telescope on Mauna Kea in Hawaii, and the SOAR Telescope at Cerro 
Pachon in Chile, with the goal of computing the object's phase function and 
nucleus size. We find an absolute magnitude of H R = 19.71 ±0.05 mag and slope 
parameter of Gr = —0.08 ± 0.05 for the inactive nucleus, corresponding to an 
effective nucleus radius of r e = 0.30 ± 0.02 km, assuming an i?-band albedo of 
Pr = 0.05. We also revisit observations reported for 259P while it was active in 
2008 to quantify the dust mass loss and compare the object with other known 
main-belt comets. 



1. Introduction 



Main-belt comets (MBCs) orbit completely within the main asteroid belt, yet exhibit 
cometary activity indicative of the sublimation of volatile ice (Hsieh & Jewitt 2006). By 
providing evidence that ice exists in the present-day asteroid belt, they present an invaluable 
opportunity to evaluate hypotheses that icy main-belt objects may have played a significant 
role in the primordial delivery of terrestrial water (e.g., Morbidelli et al. 2000). 



1 Based on observations obtained at the Southern Astrophysical Research (SOAR) telescope, which is a 
joint project of the Ministerio da Ciencia, Tecnologia, e Inovacao (MCTI) da Republica Federativa do Brasil, 
the U.S. National Optical Astronomy Observatory (NOAO), the University of North Carolina at Chapel Hill 
(UNC), and Michigan State University (MSU). 
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Comet 259P/Garradd (hereafter, 259P) was the fourth MBC to be discovered (Garradd 
et al. 2008), and has a semimajor axis of a = 2.726 AU, an eccentricity of e = 0.342, and an 



inclination of i = 15.90°. When it was active in 2008, Jewitt et al. (2009) estimated that it 
had a dust production rate on the order of 10~ 2 kg s -1 and an upper limit CN production 
rate of Qcn < 1.4 x 10 23 s -1 , and determined that it was dynamically unstable on timescales 
of ~20-30 Myr. This dynamical result indicated that 259P is likely not native to its current 
location, and probably originated elsewhere in the asteroid belt. An origin in the outer solar 
system, however, could not be ruled out. 



Jewitt et al. (2009) computed an upper limit to the effective radius of the nucleus of 



r e = 0.7 km, but since the object remained active throughout their observational campaign, 
they were unable to ascertain the true nucleus size. In this letter, we present observations 
aimed at measuring this key physical parameter. 



2. Observations 

Observations of 259P were obtained using the 8.1 m Gemini North telescope (Programs 
GN-2011A-Q-15, GN-2011B-Q-17, and GN-2012A-Q-68) on Mauna Kea in Hawaii and the 
4.2 m Southern Astrophysical Research (SOAR) telescope (Program 2012A-0461) on Cerro 
Pachon in Chile. Gemini observations were made using the imaging mode of the Gemini 



Multi-Object Spectrograph (GMOS; image scale of / .'1454 pixel" 1 ; Hook et al. 2004) and a 



Sloan Digital Sky Survey (SDSS) r' filter. SOAR observations were made using the SOAR 



Optical Imager (SOI; / .'154 pixel 1 ; Schwarz et al. 2004) and a Kron-Cousins i?-band filter. 



Non-sidereal tracking of the object was performed for all observations. 

We obtained our observations on twelve nights between UT 2011 February 28 and 2012 
May 15 during a portion of the object's orbit when active dust emission was expected to 
be minimal or absent (Table [I] Figure [T]). Inactivity was confirmed by summing the im- 
ages taken on each night and visually inspecting the resulting composite images (Figure [2]), 
ensuring that these observations reflect the properties of the inactive nucleus. 

Standard image calibration (bias subtraction and flat-field reduction) was performed for 
all images. Flat fields were constructed from dithered images of the twilight sky. Object and 
field star photometry was performed using circular apertures for which optimum sizes were 
chosen depending on the amount of field-star trailing present due to the non-sidereal tracking 
of the object, and nightly seeing conditions. Absolute calibration of object photometry 



was then performed using field star magnitudes from the SDSS catalog (York et al. 2000). 



Transformation to the BVRI system was done for all Gemini measurements, assuming solar 
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colors for the object, allowing our phase function to be computed for i?-band. 



Results 



3.1. Phase Function Analysis 

An object's brightness variation with solar phase angle (its phase function) is depen- 
dent on the surface properties of the object, such as regolith porosity, albedo, particle size 



distributions and roughness (Helfenstein & Veverka 1989 Muinonen et al. 2002). Thus, an 



understanding of the phase functions for 259P and other objects of interest, such as other 
MBCs and other asteroids and comets, provide a means for comparing their surface prop- 
erties and discerning possible contrasts or similarities. Phase function parameters also help 
us accurately determine the expected magnitude of the object at any observing geometry, 
assuming it is inactive. 

To determine the phase function of 259P, we first normalize our measured apparent 
i?-band magnitudes, m(R, A, a), to R = A = 1 AU using 



m(l, 1, a) 



m(R,A,a) - 51ogi?A 



X 



The resulting reduced magnitude, m(l, 1, a), depends only on the phase angle and rotational 
phase of the object, where the rotational properties of 259P are presently unknown. All of 
our Gemini data were taken at unknown rotational phases and thus represent "snapshot" 
observations. In our more extensive SOAR data set which covered a longer time baseline, 
we observe a photometric range (peak to trough) of Am = 0.6 ±0.1 mag. Unfortunately, we 
were unable to derive a plausible rotational lightcurve from these data, and so the rotational 
period and amplitude of 259P remain unknown. However, we see a range in brightnesses over 
our Gemini observations that is consistent with the magnitude range from our SOAR data, 
suggesting that 259P's photometric range is in fact approximately Am ~ 0.6 mag. In our 
analysis below, we assume that by using sparse sampling, we can compensate for rotational 
effects, and include such variations as a source of uncertainty in our final phase function 
solution. 



Fitting our data to the IAU H, G phase function (Bowell et al. 1989), we find best- 
fit parameters of Hr = 19.71 ± 0.05 mag and Gr = —0.08 ± 0.05. For reference, this 
corresponds to linear phase function parameters of m^(l, 1, 0) = 20.08 ± 0.05 mag and 
(3 = 0.050 ± 0.005 mag deg -1 . For comparison, the slope parameters of MBCs 133P, 176P, 
and 238P are G 
and G = 



0.04 ± 0.05 (|Hsieh et al. 2009b), G = 0.26 ± 0.05 (Hsieh et al. 2009b), 



0.03 ± 0.10 (Hsieh et al. 2011b), respectively. The range of slope parameters 
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for kilometer-scale Themis asteroids, of which 176P and 133P are identified as members, 
is —0.23 < G < 0.60 (Hsieh et al. 2008). These parameters suggest that 259P is similar 



in composition to other MBCs, which resemble dark, B-type asteroids (Hsieh et al. 2009b 



Licandro et al. 2011), and Themis family asteroids which are mostly C-type and B-type 



asteroids (Florzak et al. 1999 Tholen & Barucci 1989 Lazarro et al. 2004). 



Using a geometric i?-band albedo of pr = 0.05 (measured for other MBCs; Hsieh et al. 
2009b), we estimate an effective nucleus radius of r e ~ 0.3 km using 



p R r 2 e = (2.24 x 10 ib ) x 10 



MA[m Q -H R ] 



(2) 



where m = —27.07 is the apparent solar i?-band magnitude (Hardorp 1980 Hartmann et 



al. 1982, 1990), and r e is in km. 



3.2. Analysis of 2008 Activity 



Using the phase function we have derived for 259P, we can estimate the amount of dust 



that was present in the 2008 observations obtained by Jewitt et al. (2009). We calculate 



reduced i?-band magnitudes from their reported photometry and compare them with nu- 
cleus magnitudes that we predict from our derived phase function to determine the ratio of 
scattering surfaces of dust, A d , to that of the nucleus A N using 



1 A 



^Q0.4(m tot (l,l,a)-m fl ) 

IQ0 A(m t ot(l,l,ot)-m R ) 



(3) 



where m tot is the calculated reduced magnitude from Jewitt et al. (2009) when the total flux 



of 259P was measured, including any dust surrounding the nucleus, and tur is the expected 
reduced magnitude according to our derived IAU phase function. Using the estimated nucleus 



radius for 259P from Section |3.1[ a total dust mass, Ma, can be estimated assuming dust 
grains with radii of a = 10 pm and bulk densities of p = 1300 kg m _1 (Hsieh et al. 2004) 
using 

'A* 



ipa -A N [^) (4) 
Dust to nucleus scattering surface ratios and dust mass values derived from the reported 



photometry from Jewitt et al. (2009) are shown in Table 2 Fitting a linear function to our 



computed dust masses over this time period, we find a dust dissipation rate from the coma 
of ~ 10~ 2 kg s _1 (Figure [3]), consistent with the rate estimated by Jewitt et al. (2009). 
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Discussion 



4.1. Nature of 259P's Activity 

The discovery of disrupted asteroids, i.e., asteroids where comet-like activity is directly 
caused by an impact event and not by sublimation of volatile ices, has raised questions 
as to whether we can assume that the activity of any of the previously known MBCs is 
sublimation-driven (e.g., Jewitt||2012 ). Hsieh et al. (2012a) considered a variety of observable 
criteria that could potentially be used to distinguish between MBCs that exhibit sublimation- 
driven activity and disrupted asteroids. These criteria include long-lived activity, constant 
or increasing coma dust mass, dust morphology, and recurrent activity. 

Considering the first two criteria, we note that 259P's activity was observed over a total 
period of only 49 days and that it faded by ~60% over that period (Figure [3]). For compar- 
ison, the impact-produced dust cloud of (596) Scheila was observed for approximately one 



month, but faded by 30% in just 8 days at the end of that period (Jewitt et al. 2011)bodll. 



By contrast, the coma mass of MBC P/2010 R2 (La Sagra) increased by ~110% over 4.5 
months ( Hsieh et al.||2012c ), while the coma mass of MBC P/2006 VW139 remained constant 
for ~30 days (Hsieh et al. 2012b). The photometric behavior of the latter two objects sug- 



gests the action of prolonged dust production, consistent with sublimation and inconsistent 
with impulsive ejection. While a decrease in activity strength by itself is not a definitive 
indicator of the source of that activity, the fact that 259P's activity was only observed as it 
was fading leaves open the possibility that it was generated by an impulsive event like an 
impact, and not prolonged sublimation. 

Numerical dust modeling was used to infer the cometary nature of the activity of MBCs 
133P/Elst-Pizarro, 176P/LINEAR, 238P/Read, P/La Sagra, and P/2006 VWi 39 (|Hsieh et al. 



2004 2009a, 2011a 2012a(c Moreno et al. 2011), and also to infer the impact-driven nature 



of the activity of P/2010 A2 (LINEAR) and (596) Scheila (Jewitt et al. 2010; Snodgrass et 



al. 2010; Ishiguro et al. 2011), but has never been performed for 259P. Given the myriad 



assumptions and simplifications used in such models, however, any results would necessarily 
need to be considered inconclusive pending confirmation via other analyses (cf. Jewitt||2012 ), 
so the lack of a modeling analysis of 259P is not a critical deficiency. 

The strongest observable criterion identified to date for establishing cometary activity is 



recurrent activity with intervening periods of inactivity (Hsieh et al. 2012a Jewitt 2012). It 



is simple to reconcile such behavior with a sublimation-driven process, where activity can be 
modulated either by periodically varying solar heating over an object's orbit (as for classical 
comets), or by seasonally varying solar illumination of an isolated active site (such as the 
site of an impact where subsurface ice has been excavated and exposed; Hsieh et al.||2004" 
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Hsieh 


2009 


Capria et al. 


2012 



by other dust ejection mechanisms except in extremely contrived scenarios (Jewitt 2012) 



In the case of 259P, we have not yet had the opportunity to test for the recurrence of 
activity as it is just now about to complete one full orbit since its discovery in 2008. The 
next opportunity to observe 259P after it passes perihelion on 2013 January 25 is in July 
2013, at which point it will be at a true anomaly of v ~ 60°. In 2008, activity was observed 
between v = 29.2° and v = 48.6°. Therefore, unfortunately, it will not be possible to compare 
activity levels from 2008 and 2013 over the same orbital arcs. Nonetheless, since we have 
now measured a precise nucleus size, we can simply search for photometric excesses (beyond 
the expected brightness of a bare nucleus) in 2013 data, even if no visible cometary features 
(e.g., a tail or coma) are detected, and immediately compute the mass of any coma that may 
be present (subject to rotational phase uncertainties). Then by extrapolating the fading 
trend observed in 2008 (Table [2] Figure [3]) to the orbital arc covered by 2013 observations, 
we will be able to ascertain whether there are significant changes in activity level between 
the two epochs. 



To date, recurrent activity has only been confirmed for 133P and 238P (Hsieh et al. 



2004, 2010, 2011b). Confirming recurrent activity for 259P would strongly indicate that 
its activity is sublimation-driven, while a lack of recurrent activity would strongly suggest 
that 259P's 2008 activity could have been generated by an impact. As such, observations 
beginning in July 2013 to confirm either the presence or absence of activity are critical for 
determining this object's true nature. 



4.2. Comparison to Other MBCs 



Effective nucleus radii have now been measured (or estimated using assumed albedos) 
for five MBCs: 133P (r e = 1.9 km; |Hsieh et al. | |2009b[ |2010| |Bauer et al.l|20l2| ), 176P 
(r e = 2.0 km; |Hsieh et al.||2009b[ |2011a| |Bauer et al.||2012|), 238P (r e = 0.4 km; |Hsieh et 

r, 



al. 2011b), P/La Sagra 



0.7 km; Hsieh et al. 2012c), and now, 259P 



0.3 km; 

this work). While it is still too early to draw firm conclusions about the size distribution 
of MBCs based on measurements of just five objects, it is nonetheless interesting to note 
that all are relatively small (a few km across, or smaller), and three of the five have effective 
nucleus radii of <1 km. 

This apparent preference for small nuclei may be due to their lower escape velocities. 
Given the low ~ 1 — 3 m s _1 ejection speeds of dust particles for some MBCs (Hsieh et al. 



2004, 2009a), where the escape velocities of the largest known MBCs are ~1 m s 1 , it is 
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possible that similar emission on larger objects may simply be unobservable because ejected 
particles never actually escape the body to produce visible activity. 

Another potential explanation for the small sizes of MBC nuclei is related to the likely 
ages of these objects. The collisional destruction timescale of an object varies with its size, 
with object diameters of 10 km, 1 km, and 0.1 km corresponding to collisional destruction 
timescales of ~1-10 Gyr, ~0.1-1 Gyr, and ~10-100 Myr, respectively (Cheng 2004 Bottke 



et al.|[2005 ). Thus small objects that still exist today are likely to have formed more recently 
(e.g., in a fragmentation or cratering event involving a larger parent body) than larger ones. 
Based on their very small nucleus sizes, 238P, 259P, and P/La Sagra could therefore be much 
younger than the age of the solar system. This implication is intriguing considering indepen- 
dent dynamical findings that 133P and P/2006 VW139 may be products of extremely recent 



(<10 Myr) asteroid family-forming fragmentation events (Nesvorny et al. 2008; Novakovic 
eT^[2012| ). 



While thermal models suggest that water ice could persist for Gyrs from hundreds to 
just a few meters under the surface of a main-belt asteroid (e.g., Schorghofer|2008" Prialnik & 
Rosenberg||2009 ) , ice could be present much closer to the surface in a much younger asteroid, 



rendering that ice far more susceptible to "activation" via collisional excavation Hsieh et 



al. (cf. 2004). Extremely young asteroids would also be able to retain near-surface ice for 



wider ranges of thermal inertias, obliquities, and latitudes than those found by |Schorghofer 
(2008). As more MBCs are discovered, it will be important to continue measuring nucleus 



sizes to improve our understanding of their distribution and determine whether this apparent 
preference for small MBC nucleus sizes is real. 

In terms of activity strength relative to nucleus size (as parameterized by A d /A N ), 
259P's level of activity is similar to that of 238P (Hsieh et al. 2011b), and ~l-2 orders of 



magnitude higher than those of 133P and 176P (Hsieh et al. 2012c). However, given 259P's 



small nucleus size, in absolute terms, the total scattering cross-section of dust in the coma, 



Ad, is actually comparable to those of 133P and 176P, as well as that of 238P (Hsieh et al. 



2012c). While there are many complicating factors, such as different escape velocities for the 



wide range of MBC nucleus sizes and the fact that the amount of dust resident in the coma 
at a given time does not necessarily correspond to the instantaneous dust production rate, 
the observation that absolute activity levels are comparable among MBCs (at least within 
an order of magnitude) suggests that they may have comparably-sized active areas. If the 



collisional activation model described by Hsieh et al. (2004) is correct, this could suggest 
that the MBCs observed to date were activated by similarly-sized impactors (~1 m; Hsieh 



2009), though many more observations of many more MBCs will be needed to confirm this 



possible trend. 
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As discussed above, we were unable to determine 259P's rotational period from our 
SOAR data, which is unfortunate because a comparison to other MBCs would have been 
useful. The only MBCs with measured rotation rates — 133P (P ro t = 3.471 hr; Hsieh et 



al. 2004) and 176P (P rot = 22.23 hr; Hsieh et al. 2011a) — both exhibit extreme rotation 



(one relatively fast and one relatively slow) compared to other asteroids (cf. Masiero et al. 



2009). Our precise determinations of the size and phase function of 259P's nucleus have 



now made it possible to make much more reliable predictions of its brightness though, which 
are of obvious benefit to future observation planning. As such, attempts to measure 259P's 
rotation period at the next appropriate opportunity (i.e., when the object is expected to 
be inactive, but also close enough to the Earth and Sun to be bright enough to obtain 
high-quality photometric data) are highly encouraged. 
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Table 1. Observation Log and Photometric Results 



UT Date 


Tel. a 


N b 


t c 

^cxp 


// 


l 


R e 




^ f 


a g 


rrif/iR, A, 


a) h 


mr>( 1, 1, a) 1 


2011 


Feb 28 


Gemini 


4 


2400 


194 


5 


3 


60 


2 


71 


8 





25.4 ±0 


1 


20.5 ±0.3 


2011 


Mar 11 


Gemini 


2 


1200 


195 


8 


3 


59 


2 


77 


10 


3 


25.4 ±0 


1 


20.4 ±0.3 


2011 


Mar 26 


Gemini 


2 


1800 


197 


6 


3 


57 


2 


90 


13 


2 


25.7 ±0 


2 


20.7 ±0.4 


2011 


Mar 31 


Gemini 


2 


1800 


198 


2 


3 


57 


2 


95 


14 





25.6 ±0 


1 


20.5 ±0.3 


2011 


Apr 01 


Gemini 


4 


3600 


198 


3 


3 


57 


2 


96 


14 


1 


26.0 ±0 


1 


20.9 ±0.3 


2012 


Jan 21 


Gemini 


3 


2700 


240 


7 


2 


89 


2 


72 


19 


9 


25.4 ±0 


1 


20.9 ±0.3 


2012 


Jan 31 


Gemini 


11 


1980 


242 


5 


2 


86 


2 


55 


19 


9 


25.5 ±0 


1 


21.2 ±0.3 


2012 


Feb 01 


Gemini 


5 


900 


242 


7 


2 


86 


2 


53 


19 


9 


25.6 ±0 


1 


21.3 ±0.3 


2012 


Apr 15 


SOAR 


22 


13200 


258 





2 


59 


1 


60 


4 


4 


23.35 ±0 


05 


20.25 ±0.05 


2012 


Apr 15 


Gemini 


6 


1800 


258 





2 


59 


1 


60 


4 


4 


23.2 ±0 


1 


20.1 ±0.1 


2012 


May 13 


Gemini 


3 


540 


264 


6 


2 


49 


1 


57 


12 


4 


23.6 ±0 


1 


20.6 ±0.3 


2012 


May 15 


Gemini 


6 


1800 


265 


1 


2 


48 


1 


58 


13 


3 


23.5 ±0 


1 


20.5 ±0.3 



a Telescope 
b Number of images 

c Total effective exposure time, in seconds 
d True anomaly, in degrees 
e Heliocentric distance, in AU 
f Geocentric distance, in AU 
g Solar phase angle (Sun-259P-Earth) 
h Observed mean R-band magnitude 

'Estimated reduced R-band magnitude (normalized to R — A = 1 AU) of object at mid-point 
of full photometric range (assumed to be 0.60 mag) of rotational lightcurve 
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Table 2. Analysis of 2008 Photometry 



UT Date 


m tot {R, A,a) a 


v h 


R c 




A d 


a e 


mtotO-, l,a) f 


m R {l, 1, 


a) g 


A d /A N h 


Mi 


2008 


Sep 26 


19.69 ±0.05 


29.2 


1.85 


1 


09 


26.5 


18.16 ±0.05 


21.33 ±0 


18 


17.5 ±4.0 


8.3 ± 1.9 


2008 


Sep 30 


19.88 ±0.03 


30.9 


1.86 


1 


13 


27.1 


18.28 ±0.03 


21.35 ±0 


18 


16.0 ±3.4 


7.6 ± 1.6 


2008 


Sep 30 


19.93 ±0.03 


30.9 


1.86 


1 


13 


27.1 


18.33 ±0.03 


21.35 ±0 


18 


15.3 ±3.2 


7.3 ± 1.5 


2008 


Oct 1 


19.86 ±0.02 


31.4 


1.86 


1 


13 


27.3 


18.24 ±0.02 


21.36 ±0 


18 


16.8 ±3.3 


8.0 ± 1.6 


2008 


Oct 2 


19.90 ±0.03 


31.8 


1.87 


1 


14 


27.4 


18.26 ±0.03 


21.37 ±0 


18 


16.5 ±3.5 


7.8 ± 1.7 


2008 


Oct 3 


20.17 ±0.06 


32.3 


1.86 


1 


15 


27.6 


18.51 ±0.06 


21.37 ±0 


18 


13.0 ±3.2 


6.2 ± 1.5 


2008 


Oct 22 


21.09 ±0.04 


40.4 


1.91 


1 


33 


29.5 


19.06 ±0.04 


21.45 ±0 


19 


8.0 ±2.0 


3.8 ±0.9 


2008 


Nov 11 


21.21 ±0.10 


48.6 


1.97 


1 


56 


30.0 


18.78 ±0.10 


21.47 ±0 


19 


10.9 ±3.4 


5.2 ± 1.6 



a Observed magnitude of 259P from |Jewitt et al.| ( [2009 1 

b True anomaly, in degrees 

c Heliocentric distance, in AU 

d Geocentric distance, in AU 

°Solar phase angle (Sun-259P-Earth), in degrees 

f Total reduced magnitude (at R = A = 1 AU) of 259P 

s Expected reduced magnitude from IAU phase law 

h Inferred ratio of scattering surface area of dust to scattering surface area of the nucleus 
'Estimated dust mass within photometric aperture, in 10 kg 
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Fig. 1. — Orbital position plot for 259P observations described in Table [T] and Table [2j The 
Sun is shown at the center as a solid dot, with orbits of Mercury, Venus, Earth, Mars, 259P, 
and Jupiter (from the center outward) shown as black lines. Green squares mark positions 
where 259P was observed to be active in 2008 by Jewitt et al. (2009) and open squares 
mark inactive positions in 2011 and 2012. The perihelion (P) and aphelion (A) positions are 
marked with crosses. Observations plotted are from (a) 2008 September 26 - October 03, 
(b) 2008 October 22, (c) 2008 Nov 11, (d) 2011 February 28 - March 11, (e) 2011 March 26 
- April 01, (f) 2012 January 21 - February 01, (g) 2012 April 15, and (h) 2012 May 13-15, 
where (a)-(c) are from Jewitt et al. (2009) and (d)-(h) are from this work. 
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Fig. 2. — Composite R-band images of 259P constructed from observations obtained on 
(a) 2011 February 28 (2400 s of effective exposure time), (b) 2011 March 11 (1200 s), (c) 
2011 March 31 (1800 s), (d) 2011 April 1 (3600 s), (e) 2012 January 21 (2700 s), (f) 2012 
January 31 (1980 s), (g) 2012 April 15 (1800 s), and (h) 2012 May 15 (1800 s). All panels 
are constructed from data obtained with Gemini North. 259P is shown at the center of each 
panel measuring 30" by 30" in size, with North up and East to the left. 
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Fig. 3. — (a) Best-fit IAU phase function (solid line) for 259P, where dotted lines show the 
range of uncertainty due to phase function parameter uncertainties, and dashed lines show 
the expected photometric range due to rotational variation, assuming Am = 0.6 mag. Open 
squares mark photometry from this work (Table [T| when 259P is inactive and green squares 
indicate photometry from Jewitt et al.| ([2009) (Table [2]), in which 259P displayed cometary 
activity, (b) Estimated coma dust masses (green squares; from Table [2]) of 259P during its 
2008 active period, and a best-fit linear mass loss function (solid line), plotted versus time 
(after | Jewitt et al."]|2009[). 



